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Abstract 

We discuss supersymmetric quantum mechanical models with periodic po¬ 
tentials. The important new feature is that it is possible for both isospectral 
potentials to support zero modes, in contrast to the standard nonperiodic case 
where either one or neither (but not both) of the isospectral pair has a zero 
mode. Thus it is possible to have supersymmetry unbroken and yet also have 
a vanishing Witten index. We present some explicit exactly soluble examples 
for which the isospectral potentials have identical band spectra, and which 
are “self-isospectral” in the sense that the potentials have identical shape, but 
are translated by one half period relative to one another. 

Supersymmetry and supersymmetry breaking are fundamental issues in theoretical parti¬ 
cle physics, and supersymmetric (SUSY) quantum mechanics provides an important testing 
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ground for both physical and computational aspects of SUSY theories [TjQ]. There are also 
many applications to the theory of solitons |3[]. Of particular interest for particle physics are 
possible mechanisms for breaking SUSY dynamically. Typically, one considers models with 
discrete spectra, and then the Witten index, which characterizes the difference between the 
number of bosonic and fermionic zero modes, may be used to indicate whether or not SUSY 
is broken fl|]. Interesting subtleties arise for potentials with continuum states f|] or with 
singularities ||. 

In this paper we consider SUSY quantum mechanics for periodic potentials (which there¬ 
fore have band spectra). The main new feature is that it is possible for the periodic isospec- 
tral bosonic and fermionic potentials to have exactly the same spectrum, including zero 
modes. This is in contrast to the usual (nonperiodic and fast decaying) case for which at 
most one potential of an isospectral pair can have a zero mode. 

Consider one dimensional SUSY quantum mechanical models on the real line. The 
bosonic and fermionic Hamiltonians H± correspond to an isospectral pair of potentials V±(x) 
defined in terms of the “superpotential” W ( x ) as 

V±{x) = W 2 (x) ±W\x) (1) 

The Hamiltonians may be factorized into products of hermitean conjugate operators as 

H + = [^- + W(x)][~ + W(x)], H- = [-^ + W(x)]{± + W(x)] (2) 

which indicates that H± are formally positive operators. Thus, their energy spectrum cannot 
go below zero.[] The factorization @ also implies that V± have (almost) the same spectrum 
because there is a one-to-one mapping between the energy eigenstates : 

4+) =tM^ +w,( cK ); ^ )= M~^ +w(x) )^ (s) 

1 In some cases where V±(x) are singular H± may have negative energy states ||. We shall exclude 
such cases from the present discussion. 
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The caveat ‘almost’ is needed above because this mapping between states does not apply to 
the “zero modes” (eigenstates with E = 0), which due to the positivity of H±, are the lowest 
possible states in the spectrum. From (^) it is easy to see that the Schrodinger equation 
[—<9 2 + V±(x)]'0e = ' > = Eip^ has zero modes 

4 ± \x) = e ± J* w (4) 

provided these functions bo^ belong to the Hilbert space. SUSY is said to be unbroken if at 
least one of the bo ± ' ) is a true zero mode. Otherwise, SUSY is said to be broken dynamically. 
In the broken SUSY case there are no zero modes and so the spectra of V± are identical [due 
to the mapping (|3|)]. 

In the “standard cases” [[[[, in which V±(x ) tend to positive asymptotic values as x —> 
±oo, this means that bcf^ must be normalizable in order to be true zero modes. But it is 
clear that in these cases, at most only one of the functions bcf*^ may be normalizable. Thus 
the spectra of the potentials V± coincide except possibly for the zero energy ground state 
level. For example, if the (well-behaved) superpotential tends to asymptotic values with 
opposite signs as x —> ±oo, then one of the zero modes in ([|) is normalizable, and SUSY 
is unbroken; superpotentials that are odd functions, W(—x) = —W(x), belong to this class. 
Dynamically broken SUSY occurs when the superpotential tends to asymptotic values with 
equal signs as x —> ±oo ; superpotentials that are even in x are of this type. 

Two simple representative examples of unbroken SUSY are: (i) W(x) = x, which gives 
the harmonic oscillator. Of the two functions bo ± ' ) = e ±3;2//2 , clearly only bo f is normalizable 
and hence a zero mode, (ii) W{x) = j tanh x (with j a positive integer), which gives the 
Poschl-Tcllcr potentials 

V± = j 2 - j(j =F 1 )sech 2 x (5) 

Once again, of the two possibilities bo ± '* = [cosh x ], only bo ^ i s normalizable. U_ has j 
discrete bound states [with energies E n = n(2j — n ) for n — 0,1, ..., j — 1] and a continuum 
beginning at E — j 2 ; on the other hand, V + has j — 1 discrete bound states [with energies 
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E n = n(2j — n) for n = 1,..., j — 1] and a continuum beginning at E — j 2 . The two spectra 
coincide manifestly, except for the zero mode. 

Now consider the superpotential W(x) to be periodic, with period L: W(x + L) = W(x). 
The potentials V±{x) in (J) are therefore also periodic with period L. From the Bloch- 
Floquet theory ||, the Hilbert space consists of quasi-periodic functions: functions that 
satisfy i])k{x + L) — (exp ikL)^k(x), where the real quantity k is the crystal momentum. 

From (f|) we have F'o = ' ) (x + L) — e ± ^ L ipQ (x) where the real constant 0^ is given by 

0l= [ X+L W(y)dy (6) 

For either one of the functions 0^ t° belong to the Hilbert space, we must identify ±0^ = 
ikL. But 4>l is real, which means that 0^ = kL = 0. Thus, the two functions 0^ either 
both belong to the Hilbert space, in which case they are strictly periodic with period L : 

(x + L) — ^\x), or (when cj) L 0 0) neither of them belongs to the Hilbert space.0 
(Note that this is the exact opposite of the situation for nonperiodic potentials where if 0o ± '* 
is a zero mode of V±, then 0 q T ' ) is not a zero mode of V0.) Thus, in the periodic case the 
spectra of V + and VL match completely. 

To summarize, we see that 

0l = [ L W(y)dy = 0 (7) 

Jo 

is a necessary condition for unbroken SUSY, and when this condition is satisfied, the bosonic 
and fermionic sectors have identical spectra, including zero modes. The Witten index then 
vanishes. 

2 This conclusion is valid provided the do not have nodes (i.e. zeros, which means that the 
other function has poles). This condition is violated for example in the case of the Scarf potential 
[?j] V(x) ~ cosec 2 (x), for which W(x) = jcotx. We shall exclude such singular potentials from 
the present discussion, but note that these deserve further study - even in the nonperiodic case, 
singular superpotentials naturally exhibit interesting properties f|. 
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It is instructive to consider some simple special classes of periodic superpotentials which 
satisfy (j7j). First, suppose the superpotential is antisymmetric on a half-period: 

W{x + = -W{x) , (8) 

Then, from (11) and ([?]) we obtain 

V ± (x +1) = V r (x) (9) 

The potentials V± are simply translations of one another by half a period, and thus are 
essentially identical in shape. Therefore, they must support exactly the same spectrum, as 
SUSY indeed tells us they do. We refer to such a pair of isospectral V± that are identical in 
shape as “self-isospectral”. A simple example of a superpotential of this type is W (x) = sinx, 
with V + (x) = sin 2 (x) + cosx = U_(x + ir). 

Second, consider periodic superpotentials that are even functions of x: 

W{-x) = W(x) (10) 

but which also satisfy the condition 4>l = 0 (by subtracting an appropriate constant, any 
even W(x) can be brought into this class). The function dW(x)/dx is odd and so (jl) implies 

V±(-x) = V?(x) (11) 

The two potentials are then simply reflections of one another. They have the same shape and 
therefore give rise to exactly the same spectrum, as we know from SUSY. Such potentials are 
also “self-isospectral”. A simple example of a superpotential of this type is W(x) = cosx, 
with V + (x) = cos 2 (x) — sinx = VL (—x). 

Third and last, consider periodic superpotentials that are odd functions of x: 

W(-x) = -W(x) (12) 

Then <p L = 0 is satished trivially. The function dW(x)/dx is even and thus V±(x) are also 
even. In this case, V±(x) are not necessarily related by simple translations or reflections. 
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They are isospectral, but may not be “self-isospectral”. As an example, the superpotential 
W{x) = sinx + sin2x gives rise to an isospectral pair which is not self-isospectral, while 
W(x) = sinx + sinSx [which also belongs to the first special class mentioned above: W(x + 
7 r) = —W{x)\ gives rise to a self-isospectral pair. 

To make these general ideas more explicit, we now present a class of exactly soluble mod¬ 
els. We illustrate this class beginning with the simplest case. Consider the superpotential 


W (x) = m 


sn{x\m)cn{x\m ) 


(13) 


f4 = 


(14) 


dn(x\m) 

Here sn(x\m), cn(x\m ) and dn(x\m) are the Jacobi elliptic functions @.[|. and the (real) 
elliptic modulus parameter m can be chosen 0 < m < 1. Given this superpotential, the 
isospectral pair ([I]) of potentials is 

2 — m + 2 (m — l)/dn 2 (x\m) 

2 — m — 2dn 2 (x[m) 

Some relevant properties of the Jacobi elliptic functions are listed here. 

1. Periodicity properties: 

sn(x + 2K(rn)\m) = — sn{x\m) 
cn(x + 2K(m)\m) = — cn(x\m) 
dn{x + 2K{m)\m) = dn{x\m) 

Here K(m ) is the “real elliptic quarter period”: K{m ) = f^ 2 dd/Vl — msin 2 6. 

2. Differentiation properties: 

-^-sn(x\m) = cn(x\m) dn(x\m) 
dx 


(15) 


cn(x\m) = — sn(x\m) dn(x\m) 


d_ 

dx 

— dn(x\m) = —msn(x\m)cn(x\m) 
dx 


(16) 


3. Quadratic relations: 


— dn 2 (x\m) + 1 — m = —mcn 2 (x\m ) = m sn 2 (x\m) — m 


(17) 
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Finally, we note that when m — 1 these relations all reduce to those for the familiar 


hyperbolic functions since 


sn(cc|l) = tanhx ; cn(x |1) = sechx ; dn(x|l) = sechx 


(18) 


Thus, when m = 1 the superpotential in (|T3| ) reduces to tanhx and the isospectral potentials 
(|14]) reduce to the j — 1 case of the example in (|5|). [The reader is urged to consider the 
m —> 1 limit at all stages of the subsequent discussion.] 

From the periodicity properties (|I5|), the superpotential W(x ) in ([13]) and the potentials 


V± in (|T4|) have period 2K(m)\ 


The zero modes 


are 


V±{x + 2 K(m)) = V±{x) 


Vf \x ) = e ^ l °9Mx\m) = [ dn ( x \ m )\ 


(19) 


=Fl 


( 20 ) 


Both ^o +) and V’o ^ have period 2 K(m), and are both good zero modes. Thus the spectra 
of V± should be identical. This can be checked explicitly because the spectrum can be 
computed exactly, since the Schrodinger equation for the potentials (14) is an example of 
the Lame equation, whose explicit solution is known in terms of elliptic functions |J. Each 
spectrum has a single bound band and a continuum, as shown in Fig. 1. 

The upper edge of the bound band has energy E\ — l — m, with Bloch wavefunctions 

sn(x\m) / dn(x\m) 
cn(x\m ) 

while the the lower edge of the continuum band has energy E 2 — 1, and Bloch wavefunctions 

cn(x\m) / dn(x\m) 
sn(x\m ) 

These band-edge properties may be verified directly using the various properties listed in 

(HO- 


= 


( 21 ) 


4>t' = 


( 22 ) 


3 Note that the function dn(x\m) has no nodes or poles on the real axis. 
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Concentrating on the potential V- [an analogous analysis holds for V+\, the Schrodinger 
equation can be written [using ([T7|)] in Lame form 


ij)" = [2msn 2 (x\m) — m — E]ip 


(23) 


This equation has two independent solutions 


_ -^U ± g) Z(a) 

n ’ ~ O(x) 


(24) 


where the parameter a is related to the energy eigenvalue E by E = dn 2 (a\m ); H(x) is the 
Jacobi eta function, @(x) the Jacobi theta function, and Z(a) the Jacobi zeta function ||. 
ft is an instructive exercise to verify that at the band edges these solutions reduce to the 


wavefunctions in ( |20] . pT[ . |22|) , and furthermore that when rn — 1 they reduce to the well-known 
bound state and continuum states for the Poschl-Teller potential VL = 1 — 2 tanh 2 x. 


Given the exact solution (|24| ) we can use Bloch’s theorem to find the exact dispersion 
relation between the energy E and the crystal momentum k\ 


ip(x + 2 K(m)) = e ik2K{m) ^(x) 


k = T: 


7r 


±iZ(dn 1 (VE\m )) 


(25) 


2 K (to) 

We plotted this dispersion relation in Fig. 2, which clearly shows the band-gap. Note that 
it is rare to have an exact solution for these band features. 

The isospectral potentials V± in ( [Up are also self-isospectral. Indeed, using the properties 


sn(x + K{m)\m) = cn{x\m) / dn{x\m) 

cn(x + K(m)\m) = — \/T — m sn(x\m) / dn(x\m) 

dn(x + K (m)\m ) = \/l — m/dn(x\m) (26) 

we see that the superpotential ( |T3| ) satisfies the condition (|]), and the two potentials are 
identical up to a displacement by half a period (see Fig. |3|): 


V+(x + K(m )) = V_(x) 


( 27 ) 













Since each potential extends indefinitely and periodically, they are indistinguishable as fas 
as their spectrum is concerned. 


This raises the question of what happens in the m 1 limit, because we know that when 
m = 1 the potentials V± are genuinely different and only V- has a zero mode. The situation 
is best illustrated by Fig. 4. Consider a single period — K(m ) < x < K(m). As m —> 1, 
K{m) —> oo, and this single period becomes our real line. On this domain, the potential V- 
becomes 1 — 2 tanh 2 x, and its bound band collapses smoothly into a single discrete bound 
level (see Fig. 1). Moreover, the Bloch wavefunctions (^) and ([21]) at the lower and upper 
edges of the bound band each tend smoothly to the normalizable wavefunction i/jq ^ = sechx 
of this single bound state. On the other hand, on this domain V + flattens out and becomes 
1, which has no bound states (only a continuum E > 1). Correspondingly, the Bloch 
wavefunctions (|20| ) and (|2l|) at the lower and upper edges of the bound band tend smoothly 
to coshx and sinhx (respectively), which are not normalizable. Displacing this picture by 
K(m ) (i.e. by half a period), the roles of D_ and V + are interchanged. 

This is just the simplest example of a general class of exactly solvable periodic potentials 
with bound bands. Indeed, it is a classic result that the spectrum of the Lame equation 


ip" = \j{j + l)msn 2 px\m) — E]ip (28) 

has j bound bands and a continuum band ||. Moreover, the exact solution [analogous to 

(|4|)] can be written in terms of elliptic functions (although for j > 2 the relation between 

the energy and the crystal momentum becomes more difficult to specify explicitly). In order 

to make the connection between these Lame equations and SUSY quantum mechanics we 

must shift the Lame potential by a constant to ensure that the lower edge of the lowest band 

has energy E — 0. For example, for j = 2 the self-isospectral pair of potentials is 

{ 4m — 2 + 2 y/m 2 — m + 1 + 6(m — l)/dn 2 (x\m) 

, _ , (29) 

4 m — 2 + 2y m 2 — m + 1 — 6dn 2 [x\m) 

The energy spectrum is shown in Fig. 5, and the band-edge Bloch wavefunctions for D_ are 

E 0 = 0 : "0o ] = m + 1 + a/ m 2 — m + 1 — 3msn 2 (x[m) 
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Ei = — 1 — m + 2 \/m 2 — m + 1 : 

E 2 = — 1 + 2 m + 2 Vm 2 — m + 1 : - 02 _ ' ) 

i ?3 = 2 — m + 2 vW 2 — m + 1 : 

E^ = 4 >/m 2 — m + 1 : , 04 _ ' > 


cn(x|m)<in(a;|m) 

sn(a:|m)dn(a:|m) 

sn(a;|m)cn(a;|m) 

m + 1 — Vm 2 — m + 1 — 3msn 2 (a:|m) 


( 30 ) 


The band-edge Bloch wavefunctions -* for V+ are obtained simply by shifting the "01 ^ 
in (Si by half a period. The superpotential for the isospectral pair (§E\) is determined by 
the zero-mode 


W = ~ log 4 " } = 
ax 


6 msn(x\ 

\m)cn(x\ 

\m)dn(x\ 

\m) 

m + 1 + \Jm 

2 — m + 1 — 3 rnsn 2 (x\ 

\m) 


(31) 


When m — 1 this reduces to lb = 2 tanhx, which is the j = 2 Pbschl-Teller model in (||). 

This procedure may be repeated for higher integer values of j in the Lame equation 
(H|), leading to a general class of self-isospectral periodic potentials. The band-edge Bloch 
wavefunctions (there are 2j +1 of them since there are j bound bands) are always polynomials 
of order j in the Jacobi elliptic functions (known as Lame functions [d|). Thus, it is a 
straightforward algebraic problem to determine the band-edge wavefunctions and energies. 
Indeed, in a beautiful paper [|TD|, Alhassid et al showed that the band-edge energies are 
simply the eigenvalues of the su( 2) operator J 2 + mT 2 , with J x and J y being the standard 
su( 2 ) generators in a ( 2 j + 1 ) dimensional matrix reprentation (see also 0). 

Another generalization of the single bound band example (p~4|) is obtained by generalizing 
the superpotential ([13]) to 

sn(x\m)cn(x\m) 


W(x) = j m- 


(32) 


dn(x\m) 

where j is a positive integer. The resulting self-isospectral potentials, V± = W 2 ± W', are 


V± = j 2 (2 — m) — j(j =F l)dn 2 (x\m) + j(j ± l)(m — l)/dn 2 (x\m) (33) 

The zero mode Bloch wavefunctions are = [dn^xlm,)]^. The potentials ( p3|) are self- 
isospectral: V + {x + K(m )) = V-(x), as can be seen from Fig. 6 . Note also that as m —> 1, 
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in the domain — K(m ) < x < K(m), V- approaches the Poschl-Teller potential V- in ([5]) 
which has j discrete bound states, while V + approaches the Poschl-Teller potential V + in (||) 
which has j — 1 discrete bound states. 

We conclude by noting that it would be interesting to extend some of these ideas to held 
theoretic examples. 


Acknowledgements This work has been supported by the DOE grant DE-FG02- 
92ER40716.00 (GD) and NSF grant PHY89-04035 (JF). 


11 


REFERENCES 


[1] E. Witten, “Dynamical Breaking of Supersymmetry”, Nucl. Phys. B 185 (1981) 513, 
“Constraints on Supersymmetry Breaking”, Nucl. Phys. B 202 (1982) 253. 

[2] For a review, with references to the vast literature on SUSY QM, see: F. Cooper, 
A. Khare and U. Sukhatme, “Supersymmetry and Quantum Mechanics”, Phys. Rep. 
251 (1995) 267. 

[3] See for example: R. Rajaraman, Solitons and Instantons (North-Holland, Amsterdam, 
1982), and C. Rebbi and G. Soliani, Solitons and Particles (World Scientific, Singapore, 
1984). 

[4] R. Akhoury and A. Corntet, “Anomalous Behavior of the Witten Index - Exactly Soluble 
Models”, Nucl. Phys. B 246 (1984) 253. 

[5] A. Jevicki and J. Rodrigues, “Singular Potentials and Supersymmetry Breaking”, Phys. 
Lett. B 146 (1984) 55; J. Casahorran and S. Nam, “Singular Superpotentials and Ex¬ 
plicit Breaking of Supersymmetry”, Int. J. Mod. Phys. A 6 (1991) 2729. 

[6] C. Kittcl, Introduction to Solid State Physics (5th Edition) (Wiley, New York, 1976); 
E. C. Titchmarsh, Eigenfunction Expansions Associated with Second-Order Differential 
Equations (Vol. II) (Clarendon Press, Oxford, 1958). 

[7] F. L. Scarf, “New Soluble Energy Band Problem”, Phys. Rev. 112 (1958) 1137. 

[8] M. Abramowitz and I. Stegun (Eds.), Handbook of Mathematical Functions , (Dover, 
1990). 

[9] E. T. Whittaker and G. N. Watson, A Course of Modern Analysis, (Cambridge U. Press, 
1980); E. L. Ince, Ordinary Differential Equations , (Dover, 1956). 

[10] Y. Alhassid, F. Giirsey and F. Iachello, “Potential Scattering, Transfer Matrix, and 
Group Theory”, Phys. Rev. Lett. 50 (1983) 873. 


12 



[11] R. S. Ward, “The Nahm Equations, Finite-Gap Potentials and Lame Functions”, J. 
Phys. A: Math. Gen. 20 (1987) 2679. 

A. V. Turbiner, “Lame Equation, sis Algebra and Isospectral Deformations”, J. Phys. 
A: Math. Gen. 22 (1989) LI. 


13 



FIGURES 


i. m j 



FIG. 1. The spectrum of the isospectral pair of potentials in (|l4|), as a function of the elliptic 
parameter m. There is a single bound band, bounded below by energy Eq = 0, and above by 
energy E\ = 1 — m. There is also a continuum band beginning at E 2 = 1. Note that when m = 1 
the bound band smoothly degenerates into a single discrete bound level of energy E$ = 0; this is 
just the zero mode of the j = 1 Poschl-Teller system in (|j). 


k 



FIG. 2. The exact dispersion relation (|25|) between energy E and crystal momentum k for the 


isospectral system (14). This plot is for m = 0.3, and we clearly see the band gap between E = 0.7 
and E = 1. The horizontal line marks the edge of the Brillouin zone, at which k = ir/(2K(m)). 
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FIG. 3. The self-isospectral potentials V- (left) and V + (right). Note that they are 

identical, except for being displaced by half a period. These plots are for elliptic parameter m = 0.7. 


V_ V_ V 



FIG. 4. The upper (lower) level shows the potential VL (V+) for various values of the elliptic 
modulus parameter m = 0.5, m = 0.9, and m = 0.999, going from left to right. Notice that V+ is 
identical to VL, but shifted by half a period. Also note that as m approaches 1, within the central 
period, — K[m ) < x < K(m), V- approaches a Poschl-Teller potential with a single binding well, 
while V+ flattens out to a constant. The situation is interchanged if we displace the picture by 
K(m), which is half a period of the potentials. 
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FIG. 5. The spectrum of the isospectral pair of potentials in (p9|), as a function of the elliptic 
parameter m. There are two bound bands, and a continuum band beginning at E = 4\Jm 2 — m + 1. 
Note that when m = 1 the bound bands smoothly degenerate into two discrete bound levels of 
energy E = 0 and E = 3 and a continuum threshold at E = 4; these are the two bound states and 
continuum threshold of the j = 2 Poschl-Teller system in (|j). 
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FIG. 6. The upper (lower) level shows the potential VL (V+) in (33) for j = 3 and for various 
values of the elliptic modulus parameter m = 0.5, rn = 0.9, and m = 0.999, going from left to right. 
Notice that V+ is identical to V-, but shifted by half a period. Also note that as m approaches 1, 
within the central period, — K{m ) < x < K(m ), V- approaches the j = 3 Poschl-Teller potential, 
while V+ approaches the j = 2 Poschl-Teller potential. The situation is interchanged if we displace 
the picture by K(m), which is half a period of the potentials. 
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